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ABSTRACT 

Context. The connection of cluster mergers with the presence of extended, diffuse radio sources in galaxy clusters is still debated. 
An interesting case is the rich, merging cluster Abell 520, containing a radio halo. A recent gravitational analysis has shown in this 
cluster the presence of a massive dark core suggested to be a possible problem for the current cold dark matter paradigm. 
Aims. We aim to obtain new insights into the internal dynamics of Abell 520 analyzing velocities and positions of member galaxies. 
Methods. Our analysis is based on redshift data for 293 galaxies in the cluster field obtained combining new redshift data for 86 
galaxies acquired at the TNG with data obtained by CNOC team and other few data from the literature. We also use new photometric 
data obtained at the INT telescope. We combine galaxy velocities and positions to select 167 cluster members around z ~ 0.201. We 
analyze the cluster structure using the weighted gap analysis, the KMM method, the Dressler-Shectman statistics and the analysis of 
the velocity dispersion profiles. We compare our results with those from X-ray, radio and gravitational lensing analyses. 
Results. We compute a global line-of-sight (LOS) velocity dispersion of galaxies, o-y = 1066^^' km s"'. We detect the presence of 
a high velocity group (HVG) with a rest-frame relative LOS velocity of V|.f ~ 2000 km s"' with respect to the main system (MS). 
Using two alternative cluster models we estimate a mass range M(< 1 /i^^Mpc) = (4.0 - 9.6)xl0''' Mq. We also find that the MS 
shows evidence of subclumps along two preferred directions. The main, complex structure NGI+NG2 (with a velocity comparable 
to that of the MS) and the S'W structure (at v^f — l-l 100 km s"') define the NE-SW direction, the same of the merger suggested by 

X-ray and radio data. The £ and "W structures (at V|f 1 150 and Vrf 300 km s"') define the E-W direction. Moreover, we find 

no dynamical trace of an important structure around the lensing dark core. Rather, the HVG and a minor MS group, having different 
velocities, are roughly centered in the same position of the lensing dark core, i.e. are somewhat aligned with the LOS. 
Conclusions. We find that Abell 520 is definitely a very complex system. Our results suggest that we are looking at a cluster forming 
at the crossing of three filaments of the large scale structure. The filament aligned with the LOS and projected onto the center of the 
forming cluster might explain the apparent massive dark core shown by gravitational lensing analysis. 
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1. Introduction 

Clusters of galaxies are by now recognized to be not simple re- 
laxed structures, but rather they are evolving via merging pro- 
cesses in a hierarchical fashion from poor groups to rich clusters. 
Much progress has been made in recent years in the observations 
of the signatures of merging processes (see Feretti et al. 120021 
for a general review). A recent aspect of these investigations is 
the possible connection of cluster mergers with the presence of 
extended, diffuse radio sources: halos and relics. Cluster merg- 
ers have been suggested to provide the large amount of energy 
necessary for electron reacceleration and magnetic field amplifi- 
cation (Feretti [19991 Feretti 2002; Sarazin 2002). However, the 
question is still debated since the diffuse radio sources are quite 
uncommon and only recently we can study these phenomena 
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on the basis of a sufficient statistics (few dozen clusters up to 
z ~ 0.3, e.g., Giovannini et al. 119991 see also Giovannini & 
Feretti l2U02l Feretti ,2005j . 

Growing evidence of the connection between diffuse ra- 
dio emission and cluster merging is based on X-ray data (e.g., 
Bohringer & Schuecker 2002; Buote 2002). Studies based on 
a large number of clusters have found a significant relation be- 
tween the radio and the X-ray surface brightness (Govoni et al. 
1200 lal 1200 Ibl and connections between the presence of radio- 
halos/relics and iiTegular and bimodal X-ray surface brightness 
distribution (Schuecker et al. .2001) . 

Optical data are a powerful way to investigate the presence 
and the dynamics of cluster mergers (e.g., Girardi & Biviano 
2002), too. The spatial and kinematical analysis of member 
galaxies allow us to detect and measure the amount of sub- 
structure, to identify and analyze possible pre-merging clumps 
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or merger remnants. This optical information is really comple- 
mentary to X-ray information since galaxies and intra-cluster 
medium react on different time scales during a merger (see, e.g., 
numerical simulations by Roettiger et al. |1997l l. In this context 
we are conducting an intensive observational and data analysis 
program to study the internal dynamics of radio clusters by using 
member galaxies. Our program concerns both massive clusters, 
where diffuse radio emissions are more frequently found (e.g., 
Barrena et al. I2007bl and ref s. therein), and low-mass galaxy 
systems (Boschin et al.l200©. 

During our observational program we have conducted an in- 
tensive study of the massive cluster Abell 520 (hereafter A520). 
This cluster shows a radio halo, discovered by Giovannini et al. 
d 19991 1, having a low surface brightness with a clumpy structure 
slightly elongated in the NE-SW direction (Govoni et al. 2001bl 
[2004 , see Fig. [TTi. 

A520, al so kno wn as MS 0451H-02 in the EMSS catalog 
(Gioia et al. '1990), is a fairly rich. X-ray luminous, and hot 
cluster, with a galaxy population characterized by a high ve- 
locity dispersion: Abell richness class = 1 (Abell et al. 1989), 
Lx(0.1-2.4keV)=14.20xlO'''* h'^ erg s"' (Ebelin g et al .1996) ; 
Tx = 7.1 + 0.7 keV (Chandra da ta, Go voni et al. l2004] l: cr^ = 
(988 + 76) km s^^ (Carlberg et al.T996). 

First hints about the young dynamical status of this cluster 
came from both X-ray and optical data (Le Fevre et al. 119941 
Gioia & Luppino ' 1 994 and refs. therein). The complexity of its 
structure was confirmed by analyses of ROS AT and Chandra X- 
ray data (Govoni et al. 2001b ; 2004). In particular, new unprece- 
dent insights were recovered from deep Chandra observations by 
Markevitch et al. (I2005I I. They revealed a prominent bow shock 
indicating a cluster merger where a SW irregular structure con- 
sists of dense, cool pieces of a cluster core that has been bro- 
ken up by ram pressure as it flew in from the NE direction (see 
Fig. [TJ. The overall structure of the radio halo seems connected 
with the cluster merger and may even suggest two distinct com- 
ponents, a mushroom with a stem and a cap, where the main 
stem component goes across the cluster along the NE-SW di- 
rection and the cap ends at the bow shock (Govoni et al. 1200 Ibl 
Markevitch et al. l2005T l. 

The complex structure of A520 was also confirmed by grav- 
itational lensing analysis of Dahle et al. (2002), Mahdavi et al. 
( 120071 hereafter M07) and Okabe & Umetsu (2008). Okabe & 
Umetsu ( 120081 based on Subaru data) found a general good 
agreement between mass and galaxy luminosity distribution. 
However, the detailed study of M07 based on the same Subaru 
data and additional CFHT data pointed out a less clear situation. 
M07 found four very significant peaks in the lensing mass distri- 
bution. Among these, peaks No. 1, 2 and 4 correspond to peaks 
in the galaxy distribution and give usual values for the mass-to- 
light ratio. Peak No. 3 corresponds to the central X-ray emission 
peak, but is largely devoid of galaxies. This peak is character- 
ized by a very large mass-to-light value; thus to be referred as a 
"massive dark core". A region characterized by a somewhat low 
mass-to-light ratio exists, too (less significant peak No. 5). This 
displacement between galaxy and mass (i.e. dark matter, for the 
most part) remains very puzzling. In fact, galaxies and cold dark 
matter (CDM), being both treated as collisionless components, 
are expected to have similar behavior during a cluster merger 
If confirmed by better observations, this situation would be dif- 
ficult to explain within the widely accepted CDM paradigm of 
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cosmological structure formation (see M07 for further discus- 
sions). 

As for the analysis of the internal dynamics based on mem- 
ber galaxies, Proust et al. (I2000l l found some evidence of sub- 
structure using a sample of 21 galaxies, while the large data 
sample constructed by the Canadian Network for Observational 
Cosmology (hereafter CNOC) team (Carlberg et al. 1996; Yee et 
al. 1996) is still not exploited a part from few individual galax- 
ies in M07. Recently, we have carried out spectroscopic observa- 
tions at the TNG telescope giving new redshift data for 86 galax- 
ies in the field of A520, as well as photometric observations at 
the INT telescope. Our present analysis is based on these optical 
data as well as on the large data sample obtained by CNOC. 

This paper is organized as follows. We present our new op- 
tical data in Sect. 2 and the complete redshift catalog with the 
addition of CNOC and a few other data in Sect. 3. We present 
our results about global properties and substructure in Sect. 4. 
We furtherly analyze and discuss the dynamical status of A520 
in Sect. 5. We draw our conclusions in Sect. 6. 

Unless otherwise stated, we give eiTors at the 68% confi- 
dence level (hereafter c.l.). Throughout this paper, we use Hq = 
70 hjo km s"' Mpc"' in a flat cosmology with Qo = 0.3 and 
Qa = 0.7. In the adopted cosmology, 1' corresponds to ~ 199 
h^Q kpc at the cluster redshift. 

2. New optical data 

2.1. Spectroscopic data 

Multi-object spectroscopic observations of A520 were car- 
ried out at the TNG telescope in December 2006. We used 
DOLORES/MOS with the LR-B Grism 1, yielding a disper- 
sion of 187 A/mm, and the Loral CCD of 2048 x 2048 pixels 
(pixel size of 15 fim). This combination of grating and detec- 
tor results in a dispersion of 2.8 A/pix. We observed three MOS 
masks for a total of 102 slits. We acquired three exposures of 
1800 s for each mask. Wavelength calibration was performed us- 
ing Helium-Argon lamps. Reduction of spectroscopic data was 
carried out with the IRAF0 package. 

Radial velocities were determined using the cross- 
correlation technique (Tonry & Davis 1979) implemented in the 
RVSAO package (developed at the Smithsonian Astrophysical 
Observatory Telescope Data Center). Each spectrum was coiTe- 
lated against six templates for a variety of galaxy spectral types: 
E, SO, Sa, Sb, Sc, Ir (Kennicutt 1992). The template producing 
the highest value of fi, i.e., the parameter given by RVSAO and 
related to the signal-to-noise of the correlation peak, was cho- 
sen. Moreover, all spectra and their best correlation functions 
were examined visually to verify the redshift determination. The 
median value of H of our successfully measured galaxy redshifts 
is ~ 8. In nine cases (IDs 82, 86, 87, 143, 147 (QSO), 203, 
229, 242 and 252; see Table [B we took the EMSAO redshift 
as a reliable estimate of the redshift. Our spectroscopic survey 
in the field of A520 consists of 86 spectra with a median nom- 
inal error on cz of 60 km s The nominal eiTors as given by 
the cross-correlation are known to be smaller than the true er- 
rors (e .g., M alumuth et al. |1992| B ardelli et al. 1994; ElHngson 
& Yee 1 19941 Quintana et al. I2000I I. Double redshift determina- 
tions for the same galaxy allowed us to estimate real intrinsic 

- IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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Fig. 1. Multiwavelength picture of A520 (North is at the top and East to the left). A smoothed Chandra 0.5-2 keV image (orange 
and yellow colors) of the central region of A520 (courtesy of M. Markevitch - Markevitch et al. 2005 , X-ray point sources are 
removed) is superimposed to a r'-band image taken with the WFC camera of the INT. The contour levels of a VLA radio image at 
1 .4 GHz (courtesy of F. Govoni - Govoni et al. l2001bl l are shown, too. Main structures recovered by our analysis are highlighted (see 
Sect. |5]for more details). Label HVG indicates the center of the high velocity group having a relative LOS velocity of v,f ~ 2000 
km s"' with respect to the main system (MS). Blue circles and numbers highlight the positions of the five peaks in the lensing 
mass distribution found by M07. The size of the circles indicate the regions where we find evidence for an individual, dynamically 
important structure of the MS. The name of each structure is indicated by the label close the corresponding M07 peak number and 
the blue small square indicates the central, luminous galaxy (i.e. galaxies IDs 204, 170, 106 and 205 for N&l, N&2, S'W and 6, 
respectively). Finally, the green square indicates a head-tail radiogalaxy. 



errors in data of the same quality taken with the same instru- 
ment (Barrena et al. |2007al |2007bl l. Here we applied a similar 
correction to our nominal errors, i.e. hereafter we assume that 
true errors are larger than nominal cross-correlation errors by a 
factor 1 .5. Thus the median error on cz is 90 km s ' . 



2.2. Photometric data 

As far as photometry is concerned, our observations were car- 
ried out with the Wide Field Camera (WFC), mounted at the 



prime focus of the 2.5m INT telescope (located at Roque de los 
Muchachos observatory. La Palma, Spain). We observed A520 
in January 2008 in photometric conditions and with a seeing of 
about 2 arcsec. 

The WFC consists of a 4 chips mosaic covering a 30x30 ar- 
cmin field of view, with only a 20% marginally vignetted area. 
We took 15 exposures of 360 s using the r-SDSS (/) filter, com- 
pleting a total of 5400 s in this band. Moreover, we developed 
a dithering pattern in order to build a master "supersky" im- 
age that was used to correct our images for fringing patterns 
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Fig. 2. INT r'-band image of A520 (West at the top and North to the left). Circles and boxes indicate cluster members and non- 
member galaxies, respectively (see TablelTJ. 



(GuUixson 1 19921 1. In addition, the dithering helped us to clean 
cosmic rays and avoid gaps between CCD chips. The complete 
reduction process (including flat fielding, bias subtraction and 
bad columns elimination) yielded a final co-added image where 
the variation of the sky was lower than 1 % in the whole frame. 



Another eff'ect associated with the wide field frames is the 
distortion of the field. In order to match the photometric and 
spectroscopic samples, a good astrometric solution taking into 
account these distortions is needed. Using IRAF tasks and tak- 
ing as reference the USNO Bl.O catalog we were able to find 
an accurate astrometric solution (rms~0.3 arcsec) across the full 
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frame. The photometric caUbration was performed using Landolt 
standard fields with well known r' magnitude. These fields were 
achieved during the observation. We finally identified galaxies 
in our image and measured their magnitudes with the SExtractor 
package (Bertin & Arnouts .1996 ) and AUTOMAG procedure. 
In few cases (e.g., close companion galaxies, galaxies close to 
defects of CCD), the standard SExtractor photometric procedure 
failed. In these cases we computed magnitudes by hand. This 
method consists in assuming a galaxy profile of a typical ellipti- 
cal and scale it to the maximum observed value. The integration 
of this profile give us an estimate of the magnitude. The idea 
of this method is similar to the PSF photometry, but assuming a 
galaxy profile, more appropriate in this case. 

As a final step, we estimated and corrected the galactic ex- 
tinction, Ar' =0.12, from Burstein & Heiles (1982) reddening 
maps. We estimated that our photometric sample is complete 
down to r' = 22.0 (23.0) for S/N = 5 (3) within the observed 
field. 



agreement between our and previous redshift values. Out of the 
24 galaxies measured by Proust et al., we considered only the 
13 galaxies with > 3 and one galaxy with redshift measured 
on the emission line Wa. After having applied the correction to 
their nominal redshift, we checked the compatibility with our 
TNG-hCNOC catalog using the method described above. We 
found nine galaxies in common with our catalog for which we 
obtain mean - -0.23 + 0.46 and rms - 1 .4, in agreement with 
the expected values of and 1 . We combined TNG-i-CNOC cat- 
alog and Proust et al. data using the weighted mean of the two 
redshift determinations and the corresponding error for the nine 
galaxies in common. We added another new four galaxies by 
Proust et al., two of which are very bright galaxies. 



Table 1. Velocity catalog of 293 spectroscopically measured 
galaxies in the field of A520. In Col. 1, IDs in italics indicate 
non-cluster galaxies. 



3. Construction of the galaxy catalog 

In addition to our TNG data we considered redshifts coming 
from the CNOC survey (Carlberg et al. ^996"; Yee et al. [T996b . 
A detailed description of the data reduction techniques for the 
spectroscopic data is given in Yee et al. ( 1996). We considered 
the 215 galaxies having a redshift determined via a correlation 
significance parameter > 3 as suggested by Yee et al. (1996|l, 
see also Ellingson & Yee (1994) for the description of this pa- 
rameter As for TNG data, we applied the above correction to 
nominal errors leading to a median error on cz of ~ 120 km s ' . 
This error is quite in agreement with the error analysis performed 
by the CNOC authors (ElHngson & Yee [T994b . 

Before to proceed with the merging between TNG and 
CNOC catalogs we payed particular attention to their compati- 
bility. Twelve galaxies in the CNOC catalog are in common with 
our TNG catalog. Of these, one (galaxy ID 215) can be con- 
sidered as strongly discrepant with a Acz diff'erence of ~ 1000 
kms-' [cztng = (58982+ 100) km s-i vs.czcnoc = (5994 3+81) 
km s~']. For this galaxy a redshift by Newberry et al. ( |1988l l 
also exists and it is in agreement with the TNG redshift. For the 
remaining eleven galaxies we compared the TNG and CNOC 
determinations computing the mean and the rms of the variable 

{z\ - Zi)! sjerr^ + err^, where z\ comes from TNG, and Z2 from 
CNOC. We obtained mean = 0.53 ± 0.66 and rms = 2.2, to 
be compared with the expected values of and 1. The result- 
ing mean shows that the two sets of measurements are consis- 
tent with having the same velocity zero-point according to the 
;^^-test. The high value of rms suggests that the errors are still 
underestimated. However, when rejecting another two slightly 
discrepant determinations (Acz ~ 700 km s"' for IDs 253 and 
72) we obtained mean - 0.55 + 0.46 and rms - 1 .4, in good 
agreement with the expected values of and 1 . We decided to 
take our TNG redshifts for the galaxies IDs 215, 253 and 72 
and combine TNG and CNOC data using the weighted mean of 
the two redshift determinations and the corresponding errors for 
the remaining nine common galaxies. In total, we added another 
new 203 galaxies from CNOC obtaining a merged catalog of 289 
galaxies. 

Finally, we considered the catalog of galaxies in the field 
of A520 published by Proust et al. (2000', their Table 1). For 
three galaxies (the 2nd, 16th, and 19th) Proust et al. list only 
redshifts coming from old previous literature data. These galax- 
ies are akeady present in our TNG catalog and we verified the 
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1 


04 52 48.06, 


+03 00 43.2 


19.62 


94953 


99 


c 


2 


04 52 50.27, 


+02 57 24.7 


20.13 


68383 


153 


c 


3 


04 52 50.93, 


+02 54 37.9 


20.87 


60186 


140 


c 


4 


04 52 51.30, 


+02 58 01.7 


20.21 


98922 


104 


c 


5 


04 52 51.91, 


+02 56 44.9 


19.50 


68173 


112 


c 


6 


04 52 51.94, 


+02 59 04.0 


19.75 


108794 


225 


c 


7 


04 52 52.60, 


+02 58 21.2 


19.68 


56568 


90 


c 


8 


04 52 52.68, 


+02 54 48.6 


19.89 


60738 


117 


c 


9 


04 52 54.83, 


+02 58 54.2 


21.32 


103989 


225 


c 


10 


04 52 57.91, 


+02 54 57.3 


18.25 


60918 


108 


c 


11 


04 52 58.01, 


+02 59 16.8 


19.66 


54316 


126 


c 


12 


04 52 58.76, 


+03 00 02.7 


20.83 


90780 


94 


c 


13 


04 52 59.62, 


+02 57 34.1 


19.89 


88858 


112 


c 


14 


04 53 00.24, 


+02 57 53.3 


20.92 


99564 


90 


c 


15 


04 53 04.06, 


+02 56 09.4 


20.61 


59293 


130 


c 


16 


04 53 05.59, 


+02 56 04.1 


19.99 


64479 


126 


c 


17 


04 53 05.61, 


+03 00 35.1 


20.16 


60606 


99 


c 


18 


04 53 07.50, 


+02 57 43.4 


20.61 


116442 


225 


c 


19 


04 53 10.30, 


+03 00 12.7 


20.98 


91469 


135 


c 


20 


04 53 12.90, 


+02 55 43.5 


20.84 


77805 


135 


c 


21 


04 53 13.68, 


+02 55 05.3 


20.58 


59383 


117 


c 


22 


04 53 14.19, 


+02 55 52.0 


19.39 


61964 


108 


c 


23 


04 53 15.37, 


+02 54 19.6 


19.90 


64440 


162 


c 


24 


04 53 16.47, 


+02 56 47.3 


20.89 


98509 


126 


c 


25 


04 53 17.28, 


+02 55 15.1 


19.72 


77751 


112 


c 


26 


04 53 20.43, 


+02 58 05.0 


19.90 


60741 


130 


c 


27 


04 53 22.74, 


+02 55 59.6 


18.25 


61547 


108 


c 


28 


04 53 22.83, 


+02 55 45.4 


19.77 


62471 


99 


c 


29 


04 53 23.03, 


+02 59 18.6 


18.62 


59722 


104 


c 


30 


04 53 23.35, 


+02 57 34.8 


18.79 


62225 


112 


c 


31 


04 53 23.95, 


+02 58 33.0 


19.74 


60495 


130 


c 


32 


04 53 24.99, 


+02 57 46.2 


19.64 


78258 


126 


c 


33 


04 53 25.04, 


+03 00 25.7 


19.48 


60057 


126 


c 


34 


04 53 26.03, 


+02 56 36.7 


20.21 


98524 


130 


c 


35 


04 53 26.11, 


+02 57 46.7 


18.70 


99015 


122 


c 



In summary, our redshift catalog of A520 consists of 293 
galaxies sampling a wide, asymmetric cluster region (see Fig. [2]) 
and having a median error on cz of 112 km s"'. Table [T] lists 
the velocity catalog: identification number of each galaxy, ID 
(Col. 1); right ascension and declination, a and 6 (J2000, Col. 2); 
r' magnitudes (Col. 3); heliocentric radial velocities, v - czq 
(Col. 4) with errors, Av (Col. 5); redshift source (Col. 6; T:TNG, 
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36 


04 


53 


27 


14, 


+02 


57 


38 


8 


21 


22 


139946 


225 


C 


37 


04 


53 


29 


33, 


+02 


56 


58 


9 


19 


22 


60729 


99 


C 


38 


04 


53 


29 


61, 


+03 


00 


31 


8 


19 


66 


58588 


86 


c 


39 


04 


53 


31 


36, 


+02 


55 


09 


9 


20 


08 


79496 


130 


c 


40 


04 


53 


31 


79, 


+02 


58 


32 


1 


20 


76 


69786 


126 


c 


41 


04 


53 


32 


65, 


+02 


55 


53 


7 


18 


14 


36736 


112 


c 


42 


04 


53 


33 


63, 


+02 


54 


57 


9 


19 


61 


63283 


140 


c 


43 


04 


53 


34 


62, 


+02 


56 


32 


4 


19 


58 


59880 


99 


c 


44 


04 


53 


35 


76, 


+02 


58 


31 


5 


17 


29 


59488 


122 


c 


45 


04 


53 


36 


05, 


+02 


55 


03 


4 


20 


86 


59353 


126 


c 


46 


04 


53 


36 


54, 


+03 


00 


01 


6 


21 


03 


59401 


144 


c 


47 


04 


53 


36 


76, 


+02 


56 


37 


3 


18 


85 


64659 


117 


c 


48 


04 


53 


36 


99, 


+02 


57 


47 


3 


19 


84 


59320 


99 


c 


49 


04 


53 


37 


01, 


+02 


54 


48 


7 


20 


90 


79073 


104 


c 


50 


04 


53 


38 


36, 


+02 


57 


31 


6 


19 


88 


59272 


130 


c 


51 


04 


53 


39 


06, 


+02 


57 


10 


3 


18 


99 


59707 


112 


c 


52 


04 


53 


41 


08, 


+02 


58 


09 


1 


20 


85 


62576 


94 


c 


53 


04 


53 


41 


56, 


+02 


55 


23 


4 


18 


18 


62962 


108 


c 


54 


04 


53 


41 


88, 


+02 


57 


29 


9 


19 


84 


59832 


135 


c 


55 


04 


53 


41 


97, 


+02 


59 


00 


3 


19 


26 


60168 


130 


c 


56 


04 


53 


42 


44, 


+02 


55 


09 


4 


19 


68 


58241 


126 


c 


57 


04 


53 


42 


56, 


+02 


57 


33 


9 


18 


73 


64665 


122 


c 


58 


04 


53 


42 


86, 


+02 


59 


58 


9 


18 


55 


44786 


126 


c 


59 


04 


53 


42 


89, 


+02 


54 


19 


3 


20 


37 


58729 


117 


c 


60 


04 


53 


43 


19, 


+03 


00 


04 


2 


20 


18 


98703 


158 


c 


6i 


04 


53 


43 


32, 


+02 


58 


59 


9 


18 


67 


65472 


94 


c 


62 


04 


53 


43 


78, 


+02 


59 


22 


2 


20 


96 


79478 


94 


c 


63 


04 


53 


43 


81, 


+02 


56 


12 


9 


18 


65 


65247 


130 


c 


64 


04 


53 


46 


26, 


+02 


56 


45 


8 


18 


68 


64569 


108 


c 


65 


04 


53 


46 


77, 


+02 


54 


19 


9 


20 


07 


61199 


158 


c 


66 


04 


53 


47 


93, 


+02 


57 


30 


2 


20 


05 


59578 


117 


c 


67 


04 


53 


49 


09, 


+02 


55 


50 


5 


18 


97 


75509 


108 


c 


68 


04 


53 


49 


45, 


+02 


56 


50 


2 


20 


27 


65478 


126 


c 


69 


04 


53 


50 


05, 


+02 


55 


02 


7 


18 


64 


59467 


135 


c 


70 


04 


53 


50 


16, 


+02 


57 


10 


2 


19 


60 


61436 


122 


c 



iU 






Q-,(5(J2000) 






r 




V 

( Km 


ZiV 

s ) 


Source 


71 


04 


53 


52 


22, 


+02 59 


32 


5 


21 


78 


65364 


117 


C 


72 


04 


53 


52 


44, 


+02 54 


14 


7 


19 


91 


65804 


58 


T 


73 


04 


53 


52 


53, 


+02 54 


32 


7 


19 


57 


62582 


130 


C 


74 


04 


53 


53 


31, 


+02 55 


04 


5 


20 


59 


115810 


225 


C 


75 


04 


53 


53 


46, 


+02 52 


12 


1 


18 


86 


67466 


102 


T 


76 


04 


53 


54 


33, 


+02 56 


39 


2 


19 


29 


61182 


126 


C 


77 


04 


53 


55 


31, 


+02 55 


56 


5 


18 


26 


60198 


122 


C 


78 


04 


53 


55 


74, 


+02 57 


46 


6 


18 


83 


59254 


126 


C 


79 


04 


53 


55 


99, 


+02 48 


45 


7 


19 


21 


61134 


60 


T 


80 


04 


53 


57 


31, 


+02 52 


54 


4 


19 


44 


62808 


135 


T 


81 


04 


53 


58 


20, 


+03 00 


37 


8 


19 


49 


59476 


140 


C 


82 


04 


53 


58 


33, 


+02 47 


18 


1 


19 


47 


115868 


64 


T 


83 


04 


53 


58 


35, 


+02 58 


30 


5 


19 


54 


60492 


748 


C 


84 


04 


53 


58 


54, 


+02 49 


45 





19 


10 


65115 


87 


T 


85 


04 


53 


58 


83, 


+02 59 


33 


4 


19 


34 


60357 


99 


C 


86 


04 


53 


59 


04, 


+02 50 


57 


3 


19 


67 


64947 


27 


T 


87 


04 


53 


59 


29, 


+02 53 


07 


2 


20 


41 


76130 


21 


T 


88 


04 


53 


59 


36, 


+02 51 


17 


1 


19 


45 


76190 


63 


T 


89 


04 


53 


59 


59, 


+02 56 


36 


5 


20 


04 


62429 


148 


C 


90 


04 


53 


59 


99, 


+02 59 


45 


1 


18 


40 


49274 


117 


C 


91 


04 


54 


00 


34, 


+03 03 


55 


9 


19 


39 


70008 


80 


T 


92 


04 


54 


00 


56, 


+02 53 


32 


9 


19 


58 


61319 


144 


C 


93 


04 


54 


00 


73, 


+02 48 


19 


2 


19 


66 


64849 


88 


T 


94 


04 


54 


00 


88^ 


+02 59 


16 


5 


20 


62 


77280 


135 


C 


95 


04 


54 


01 


15, 


+02 57 


45 


6 


17 


35 


62154 


51 


C + P 


96 


04 


54 


01 


76, 


+02 55 


47 


4 


19 


17 


58327 


90 


C 


97 


04 


54 


02 


37, 


+03 01 


58 


2 


19 


34 


59542 


68 


T 


98 


04 


54 


02 


49, 


+02 59 


31 





19 


25 


49364 


320 


C 


99 


04 


54 


02 


71, 


+02 50 


36 


2 


19 


61 


58860 


68 


T 


100 


04 


54 


02 


88, 


+02 52 


22 


6 


18 


63 


60818 


76 


T 


101 


04 


54 


02 


90, 


+02 51 


06 


7 


19 


61 


60461 


116 


T 


102 


04 


54 


03 


05, 


+02 49 


34 


7 


20 


24 


58347 


117 


T 


103 


04 


54 


03 


26, 


+03 04 


40 


5 


20 


83 


152875 


112 


T 


104 


04 


54 


03 


36, 


+02 55 


40 


3 


19 


74 


70034 


144 


C 


105 


04 


54 


03 


45, 


+02 59 


30 


6 


18 


79 


60516 


122 


C 



C:CNOC and P:Proust et al.). We list r' magnitudes for 291 out 
of 293 galaxies having redshifts. The exceptions are a galaxy 
just outside the western border of the imaging field and a huge 
foreground spiral galaxy. We have redshifts for galaxies down 
to r' ~21.5 mag, but we are 40% complete down to r' = 19 
mag within 3 arcmin from R.A.=04*'54'"14% Dec.=+02°57'00" 
(J2000.0). The completeness of the spectroscopic sample de- 
creases in the outskirts of the cluster. 

Figure[3]shows the contribute of TNG data added to previous 
spectroscopic information. 

A520 does not exhibit the presence of a clear dominant 
galaxy and in fact it is classified as Bautz-Morgan class III 
(Abell et al. 1989). In particular, our sample lists nine luminous 
galaxies in a range of one mag from the most luminous one: 
IDs 204, 106, 44, 170, 95, 205, 264, 160 and 137. These galax- 
ies are generally sparse in the field. A few of these galaxies are 
close to the lensing mass peaks pointed out by M07, i.e. ID 204 
is close to peak No. 1; the galaxy couple composed by IDs 160 
and 170 is close to peak No. 2; ID 106 is close to peak No. 4; 
ID 205 is close to peak No. 5. 

Govoni et al. 1* 200 Ibl ) pointed out the presence of several dis- 
crete radio sources in the field of A520. In particular, there are 
two head-tail radio sources (0454+0255 A and 0454+0255B; see 
also Cooray et al. il998 ) located on the eastern side with the tails 



oriented toward the same direction, opposite to the cluster center. 
Our catalog lists the redshift for the northern one, 0454+0255B 
(ID 184), which is classified as a cluster member. Cooray et al. 
(119981 ) also list a third radio source (0454+0257) which, again, is 
classified as a cluster member (ID 95). From a visual inspection 
of the Chandra image studied by Markevitch et al. (120051 Obs.Id 
42 1 5) we also note that 0454+0257 is an evident pointlike X-ray 
source in the field of A520. 

4. Analysis and results 

4.1. Member selection 

To select cluster members out of 293 galaxies having redshifts, 
we follow a two steps procedure. First, we perform the adaptive- 
kernel method (hereafter DEDICA, Pisan i 1993 and 1996; see 
also Fadda et al. 119961 Girardi et al. 11996 ; Girardi & Mezzetti 
.200 Ij) . We search for significant peaks in the velocity distribu- 
tion at >99% C.I.. This procedure detects A520 as an asymmet- 
ric one-peak structure at z ~ 0.201 populated by 223 galaxies 
considered as candidate cluster members (see Fig.|4]i. 

All the galaxies assigned to the A520 peak are analyzed in 
the second step which uses the combination of position and ve- 
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ID 


Q',(5(J2000) 


r' 


V 


Av 


Source 










(1cm 


s-') 




106 


04 54 03.82, 


+02 53 32.4 


17.10 


61277 


114 


T 


107 


04 54 03 96 


+02 53 40.7 


17.98 


64646 


70 


T + C 


108 


04 54 04.18, 


+03 02 48.4 


20.84 


55511 


104 


T 


109 


04 54 04.30, 


+02 49 00.8 


20.01 


59783 


141 


T 


110 


04 54 04 54 


+02 52 43.4 


19.08 


60827 


98 


T 


111 


04 54 04.59, 


+02 56 54.2 


20.12 


61145 


112 


P 


112 


04 54 04 67 


+02 56 04 


20.18 


60327 


148 


C 


113 


04 54 05 13 


+02 47 09.1 


19.26 


65020 


80 


T 


114 


04 54 05 14 


+02 56 22.2 


20.02 


59536 


130 


C 


115 


04 54 05.38, 


+03 04 31.0 


20.26 


70903 


92 


T 


116 


04 54 05 44 


+02 59 15 5 


19.76 


66554 


117 


C 


117 


04 54 05.92, 


+02 55 54.2 


19.92 


59916 


72 


C + P 


118 


04 54 05.92, 


+02 55 46.0 


20.15 


58669 


108 


C 


119 


04 54 05.93, 


+02 53 37.3 


19.82 


59458 


153 


C 


120 


04 54 06 02 


+02 57 56 7 


21.22 


62839 


104 


C 


121 


04 54 06.12, 


+02 58 46.1 


20.20 


62018 


450 


C 


122 


04 54 06 32 


+03 03 27 9 

\\J^ \j.y 1 ty 


19.67 


59639 


128 


T 


123 


04 54 06.64, 


+02 59 06.2 


20.09 


98850 


122 


C 


124 


04 54 06 77 


+02 53 55 4 


19.59 


61473 


71 


T + C 


125 


04 54 06.78, 


+02 57 41.6 


20.39 


66990 


54 


C + P 


126 


04 54 06 82 


+03 03 56 

\ \J-y \J-y my\jt\j 


19.38 


70420 


64 


T 


127 


04 54 07.47, 


+02 51 44.6 


20.70 


59358 


93 


T 


128 


04 54 07 62 


+03 00 59 2 


19.31 


59490 


46 


T + C 


129 


04 54 07.74, 


+02 56 00.4 


19.78 


60861 


104 


C 


130 


04 54 07.83, 


+02 57 02.7 


19.82 


61140 


122 


C 


131 


04 54 08.47, 


+02 59 08.4 


20.03 


60060 


117 


C 


132 


04 54 08.63, 


+02 56 36.1 


19.19 


62695 


94 


C 


133 


04 54 08.78, 


+03 01 21.4 


19.10 


60931 


72 


T 


134 


04 54 08.87, 


+02 53 49.5 


18.69 


57872 


117 


C 


135 


04 54 08.90, 


+02 53 21.5 


19.11 


59219 


66 


T 


136 


04 54 09.03, 


+02 52 00.8 


18.63 


63622 


66 


T 


137 


04 54 09.06, 


+02 59 48.8 


17.93 


60075 


76 


C 


138 


04 54 09.37, 


+02 55 15.8 


18.97 


61125 


117 


C 


139 


04 54 09.41, 


+02 50 22.7 


19.66 


60440 


104 


T 


140 


04 54 09.41, 


+02 51 32.1 


19.38 


61024 


48 


T 



iU 


Q-,(5(J2000) 


r' 


V 




Source 










( Km 


s ) 




141 


04 54 09.42, 


+02 56 26.3 


19.89 


67036 


94 


C 


142 


04 54 09.55, 


+02 55 40.5 


19.20 


58879 


112 


C 


143 


04 54 09.60, 


+03 02 22.9 


20.83 


181270 


51 


T 


144 


04 54 10.10, 


+02 55 42.2 


19.69 


60600 


117 


C 


145 


04 54 10.31, 


+02 54 38.9 


19.41 


60963 


51 


T + C 


146 


04 54 10.41, 


+02 56 09.9 


20.02 


58280 


148 


C 


147 


04 54 10.52, 


+02 47 39.2 


19.83 


704502 


100 


T 


148 


04 54 10.69, 


+03 02 20.0 


20.11 


152100 


96 


T 


149 


04 54 11.48, 


+02 55 25.8 


20.26 


57941 


108 


C 


150 


04 54 11.69, 


+02 59 13.1 


19.98 


62291 


135 


C 


151 


04 54 11.79, 


+02 48 10.7 


18.28 


58729 


56 


T 


152 


04 54 11.80, 


+02 52 11.4 


19.60 


60376 


54 


T 


153 


04 54 11.82, 


+02 50 48.1 


19.11 


59577 


182 


T 


154 


04 54 11.93, 


+02 58 07.8 


18.23 


62292 


47 


C + P 


155 


04 54 12.08, 


+02 56 36.5 


20.77 


60708 


148 


C 


156 


04 54 12.19, 


+02 57 50.7 


20.48 


58417 


153 


C 


157 


04 54 12.31, 


+03 02 47.8 


20.89 


85271 


188 


T 


158 


04 54 12.77, 


+02 49 56.5 


19.66 


82005 


78 


T 


159 


04 54 13.04, 


+02 56 33.2 


19.36 


60084 


144 


C 


160 


04 54 13.14, 


+02 57 33.8 


17.70 


60115 


24 


P 


161 


04 54 13.16, 


+02 58 36.6 


20.25 


62495 


135 


C 


162 


04 54 13.34, 


+03 02 08.8 


20.65 


99114 


106 


T 


163 


04 54 13.35, 


+02 51 58.1 


20.05 


60776 


96 


T 


164 


04 54 13.50, 


+02 48 33.7 


20.45 


75930 


180 


T 


165 


04 54 13.68, 


+02 56 10.2 


19.78 


59653 


99 


C 


166 


04 54 13.74, 


+02 53 26.7 


19.14 


60519 


72 


T 


167 


04 54 13.80, 


+02 59 19.4 


19.81 


59059 


130 


C 


168 


04 54 14.01, 


+02 55 42.5 


19.63 


59383 


122 


C 


169 


04 54 14.09, 


+03 01 05.1 


20.39 


40586 


144 


C 


170 


04 54 14.10, 


+02 57 09.9 


17.29 


59506 


69 


P 


171 


04 54 14.17, 


+03 01 10.3 


18.89 


66961 


140 


C 


111 


04 54 14.34, 


+02 58 36.5 


18.50 


60111 


104 


C 


173 


04 54 14.36, 


+02 59 16.3 


19.41 


58513 


117 


C 


174 


04 54 14.40, 


+02 56 42.2 


19.02 


60762 


99 


C 


175 


04 54 14.79, 


+03 00 49.0 


19.01 


58513 


75 


T + C 



locity information: the "shifting gapper" method by Fadda et al. 
( I1996I I. This procedure rejects galaxies that are too far in veloc- 
ity from the main body of galaxies and within a fixed bin that 
shifts along the distance from the cluster center. The procedure 
is iterated until the number of cluster members converges to a 
stable value. Following Fadda et al. (|1996) we use a gap of 1000 
km s"' - in the cluster rest-frame - and a bin of 0.6 /i^,,' Mpc, 
or large enough to include 15 galaxies. The choice of the cluster 
center is not obvious. In fact, several galaxy condensations are 
visible in the field (Gioia & Luppino 1994). Moreover, no obvi- 
ous dominant galaxy is present (see Sect. [3]) and the lensing mass 
distribution shows several peaks (e.g., M07). Thus, hereafter we 
assume the position of the peak of X-ray emission as listed by 
Ebeling et al. (1996) [R.A.=04'>54'"07 ?44, Dec.=+02°55'12 :'0 
(J2000.0)] as the cluster center. After the "shifting gapper" pro- 
cedure we obtain a sample of 167 fiducial cluster members (see 
Fig.E). 

The 2D galaxy distribution analyzed through the 2D 
DEDICA method shows only one peak [at R.A. =04^54™ 13 ?55, 
Dec.=+02°56'35 :'2 (J2000.0)]. This peak, hereafter the "opti- 
cal" cluster center, is displaced towards NE with respect to the 
X-ray peak and is close, but not coincident, to a pair of luminous 
galaxies (IDs 160 and 170). The biweight cluster center, i.e. that 
recovered by computing the biweight means (Beers et al. 1 19901 1 



of R.A. and Dec. of galaxy positions [R.A.=04h54'"12 ?62, 
Dec.=+02°55'57:'3 (J2000.0)], is roughly coincident with the 
DEDICA peak. Using these alternative cluster centers we verify 
the robustness of our member selection. 

4.2. Global klnematlcal properties 

By applying the biweight estimator to the 167 cluster mem- 
bers (Beers et al. I1990I I. we compute a mean cluster redshift 
of <z) = 0.2008+ 0.0003, i.e. (v) = (60209+82) km s-\ We 
estimate the LOS velocity dispersion, o-y, by using the biweight 
estimator and applying the cosmological correction and the stan- 
dard correction for velocity eiTors (Danese et al. I1980I I. We ob- 
tain cTy - 1066^J^j km s"', where eiTors are estimated through a 
bootstrap technique. 

To evaluate the robustness of the o-y estimate we analyze the 
velocity dispersion profile (Fig.|6]l. The integral profile smoothly 
decreases and flattens beyond ~ 0.6 hj^^ Mpc suggesting that 
a robust value of o-y is asymptotically reached in the external 
cluster regions, as found for most nearby clusters (e.g., Fadda et 
al. [19961 Girardi et al. lT996] l. 
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ID 



a,6(}2000) 



V 

(km 



Av Source 



176 


04 


54 


15 


09, 


+02 


57 


07 


8 


18 


06 


59163 


80 


C + P 


177 


04 


54 


15 


55, 


+02 


54 


58 


3 


20 


26 


18338 


270 


c 


178 


04 


54 


15 


76, 


+02 


52 


46 


9 


19 


09 


62064 


122 


T 


179 


04 


54 


15 


89, 


+03 


04 


47 


6 


20 


62 


58757 


122 


T 


180 


04 


54 


15 


95, 


+02 


58 


19 


1 


19 


43 


60267 


135 


C 


181 


04 


54 


16 


01, 


+02 


55 


20 


7 


18 


33 


60954 


104 


C 


182 


04 


54 


16 


06, 


+02 


56 


42 


8 


18 


64 


58821 


94 


C + P 


183 


04 


54 


16 


56, 


+02 


57 


26 


7 


19 


67 


58876 


126 


C 


184 


04 


54 


16 


57, 


+02 


55 


31 


8 


19 


35 


60972 


104 


C 


185 


04 


54 


16 


89, 


+02 


54 


24 


8 


20 


25 


62827 


103 


C 


186 


04 


54 


16 


94, 


+02 


48 


37 


4 


20 


03 


82173 


96 


T 


187 


04 


54 


17 


10, 


+03 


01 


49 


5 


19 


35 


60006 


57 


T 


188 


04 


54 


17 


31, 


+02 


53 


12 





19 


07 


60264 


60 


T 


189 


04 


54 


17 


33, 


+02 


56 


46 


1 


19 


60 


64056 


122 


P 


190 


04 


54 


17 


43, 


+02 


59 


24 





19 


54 


58780 


76 


C 


191 


04 


54 


17 


66, 


+02 


48 


24 


9 


19 


66 


82097 


78 


T 


192 


04 


54 


17 


90, 


+02 


55 


35 





19 


31 


60549 


117 


C 


193 


04 


54 


17 


95, 


+02 


46 


49 


6 


18 


94 


60927 


72 


T 


194 


04 


54 


18 


02, 


+02 


57 


41 


5 


20 


14 


61397 


126 


C 


195 


04 


54 


18 


18, 


+02 


59 


55 


7 


20 


09 


61334 


130 


C 


196 


04 


54 


18 


58, 


+03 


00 


36 


5 


20 


45 


57385 


90 


T 


197 


04 


54 


18 


88, 


+02 


50 


54 


4 


16 


93 


18754 


100 


T 


198 


04 


54 


19 


00, 


+02 


56 


17 


2 


18 


74 


62504 


135 


C 


199 


04 


54 


19 


05, 


+02 


56 


13 


8 


20 


94 


91532 


130 


C 


200 


04 


54 


19 


16, 


+02 


58 


26 


5 


19 


10 


60552 


117 


C 


201 


04 


54 


19 


28, 


+03 


01 


09 


9 


20 


64 


58353 


122 


T + C 


202 


04 


54 


19 


31, 


+02 


51 


47 


5 


19 


71 


58948 


60 


T 


203 


04 


54 


19 


51, 


+02 


48 


05 


7 


19 


35 


8555 


31 


T 


204 


04 


54 


19 


91, 


+02 


57 


44 


8 


16 


93 


60315 


64 


C + P 


205 


04 


54 


19 


96, 


+02 


55 


30 


6 


17 


35 


58597 


99 


C 


206 


04 


54 20 


17, 


+02 


55 


32 


5 


19 


43 


58381 


108 


C 


207 


04 


54 20 


21, 


+02 


59 


20 


9 


20 


88 


59955 


176 


C 


208 


04 


54 20 


56, 


+03 


00 


55 


9 


19 


94 


76010 


123 


T + C 


209 


04 


54 20 


58, 


+02 


53 


37 


4 


19 


52 


59946 


90 


T + C 


210 


04 


54 20 


62, 


+02 


56 


41 


4 


19 


85 


58657 


81 


C 



ID 



Q-, (5 (12000) 



Av Source 

























(km 
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211 


04 


54 


20 


68, 


+02 


55 


29 


8 


18 


32 


58969 


104 


C 


212 


04 


54 


21 


07, 


+02 


51 


24 


9 


19 


08 


61484 


72 


T 


213 


04 


54 


21 


72, 


+02 


55 


56 





20 


03 


58244 


140 


C 


214 


04 


54 


21 


73, 


+03 


05 


11 


6 


19 


79 


61250 


74 


T 


215 


04 


54 


21 


84, 


+02 


55 


00 





19 


78 


58982 


100 


T 


216 


04 


54 


22 


02, 


+02 


57 


43 


3 


22 


27 


110602 


225 


C 


217 


04 


54 


23 


08, 


+02 


50 


18 


9 


20 


93 


165148 


88 


T 


218 


04 


54 


23 


13, 


+02 


58 


01 


2 


20 


03 


60228 


130 


C 


219 


04 


54 


23 


27, 


+02 


57 


09 









17085 


108 


C 


220 


04 


54 


23 


27, 


+02 


59 


13 





19 


83 


61202 


94 


C 


221 


04 


54 


23 


48, 


+03 


03 


16 


1 


20 


94 


78236 


238 


T 


222 


04 


54 


23 


53, 


+02 


50 


34 


5 


18 


50 


60667 


69 


T 


223 


04 


54 


23 


84, 


+02 


51 


08 


6 


20 


14 


114162 


98 


T 


224 


04 


54 


23 


99, 


+02 


56 


10 


9 


21 


16 


70466 


108 


C 


225 


04 


54 


24 


24, 


+03 


01 


11 


9 


21 


50 


79751 


117 


T 


226 


04 


54 


24 


88, 


+02 


58 


56 





19 


93 


60289 


87 


C + P 


227 


04 


54 


24 


95, 


+02 


52 


23 


5 


20 


02 


60260 


90 


T 


228 


04 


54 


25 


34, 


+02 


47 


03 


9 


19 


70 


61235 


72 


T 


229 


04 


54 


25 


37, 


+02 


49 


23 


4 


20 


62 


57854 


226 


T 


230 


04 


54 


25 


50, 


+02 


59 


38 


3 


18 


98 


60057 


104 


C 


231 


04 


54 


26 
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+02 


54 


23 


3 


21 


25 


50350 
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C 


232 


04 


54 


26 
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+03 


00 


46 


9 


19 


87 


60182 


77 
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04 


54 
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80, 


+02 


58 


21 


9 


19 


86 


78677 
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C 
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04 


54 


27 
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+03 


03 


29 


3 


20 


02 


60186 


99 


T 


235 


04 


54 


27 


76, 


+02 


55 


29 


2 


18 


63 


67171 


90 


C 


236 


04 


54 


27 


96, 


+02 


54 


18 


1 


19 


61 


61070 


93 


T 


237 


04 


54 


28 


14, 


+02 


55 


45 


7 


20 


41 


59611 


112 


C 


238 


04 


54 


28 


18, 


+02 


55 


36 


5 


18 


74 


66893 


99 


C 


239 


04 


54 


28 


63, 


+03 


04 


16 





19 


17 


60508 


62 


T 


240 


04 


54 


29 


02, 


+02 


54 


29 





19 


97 


65186 


61 


T + C 


241 


04 


54 


29 


02, 


+02 


56 


59 


5 


20 


54 


62384 


153 


C 


242 


04 


54 


29 


14, 


+02 


48 


55 


9 


19 


74 


135190 


366 


T 


243 


04 


54 


29 


53, 


+02 


58 


22 


1 


19 


34 


67096 
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C 


244 


04 


54 


29 


58, 


+02 


55 


22 





19 


66 


62063 


112 


C 


245 


04 


54 


30 


21, 


+03 


02 


25 


3 


18 


59 


66861 


98 


T 



4.3. Substructure 

4.3.1. Velocity distribution 

We analyze the velocity distribution to look for possible devia- 
tions from Gaussianity that might provide important signatures 
of complex dynamics. For the following tests the null hypothesis 
is that the velocity distribution is a single Gaussian. 

We estimate three shape estimators, i.e. the kurtosis, the 
skewness, and the scaled tail index (see, e.g.. Beers et al. |1991| l. 
According to the value of the skewness (+0.471) the velocity 
distribution is positively skewed and differs from a Gaussian at 
the 95 - 99% c.l. (see Table 2 of Bird & Beers 1993 ). Moreover, 
according the the scaled tail index the velocity distribution is 
heavily tailed and differs from a Gaussian at the 90 - 95% c.l. 
(see Table 2 of Bird & Beers [T993T l. 

Then we investigate the presence of gaps in the velocity dis- 
tribution. A weighted gap in the space of the ordered velocities 
is defined as the difference between two contiguous velocities, 
weighted by the location of these velocities with respect to the 
middle of the data. We obtain values for these gaps relative to 
their average size, precisely the midmean of the weighted-gap 
distribution. We look for normalized gaps larger than 2.25 since 
in random draws of a Gaussian distribution they arise at most in 



about 3% of the cases, independent of the sample size (Wainer 
and Schacht |1978l see also Beers et al. |1991| l. We detect a sig- 
nificant gap (at the 99.95% c.l.) which separates the main cluster 
from a group of 22 high velocity galaxies (see Fig.|5]and the first 
line of Table|2]i. For each gap Table|2]lists the number of galaxies 
for the group before the gap and that after the gap (Col. 2); the 
velocity boundaries before and after the gap (Col. 3); the size 
of the gap (Col. 4); the probability of finding such a gap in a 
Gaussian distribution (Col. 5). Hereafter we define MS the main 
system with the 145 galaxies having low velocities and HGV the 
group with the 22 galaxies having high velocities (see Table [3] 
for their main kinematical properties). 

As for the spatial distribution, there is no difference be- 
tween the galaxies of the HVG and the MS (according to the 
2D Kolmogorov-Smirnov test - hereafter 2DKS-test - Fasano 
& Franceschini 119871 see Fig. |2ll. However, when considering 
the clustercentric distances, they differ at the 97.5% c.l. accord- 
ing to the ID Kolmogorov-Smirnov test (hereafter IDKS-test, 
see e.g.. Press et al. 1992) with the HVG galaxies being, on aver- 
age, closer to the cluster (X-ray) center Accordingly, while the 
optical (biweight) center of the MS lies very close to the optical 
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246 


04 54 30.31 


+02 58 44.8 


19.18 


59877 


104 


C 


247 


04 54 30.67 


+02 54 44.4 


20.05 


66872 


117 


C 


248 


04 54 30.90 


+02 59 35.4 


19.93 


17451 


225 


c 


249 


04 54 31.01 


+02 49 04.3 


20.05 


60211 


64 


T 


250 


04 54 31.23 


+03 05 13.6 


18.97 


60001 


88 


T 


251 


04 54 31.42 


+02 57 22.9 


20.36 


74618 


180 


C 


252 


04 54 31.93 


+02 52 36.7 


20.39 


76231 


41 


T + C 


253 


04 54 32.31 


+03 03 52.1 


20.23 


59672 


147 


T 


254 


04 54 32.63 


+02 53 01.3 


20.32 


59221 


164 


T 


255 


04 54 32.68 


+02 54 48.9 


18.00 


60894 


59 


C + P 


256 


04 54 33.56 


+03 03 23.1 


19.10 


99246 


112 


C 


257 


04 54 33.78 


+02 58 51.8 


20.91 


125307 


225 


C 


258 


04 54 34.25 


+02 50 00.3 


19.46 


111247 


81 


T 


259 


04 54 35.27 


+03 01 05.3 


18.84 


59269 


104 


C 


260 


04 54 35.83 


+03 01 05.4 


20.53 


61005 


153 


C 


261 


04 54 36.49 


+02 54 33.7 


21.45 


112128 


225 


C 


262 


04 54 36.93 


+03 03 23.7 


19.63 


60798 


122 


C 


263 


04 54 37.47 


+03 02 24.4 


18.65 


60021 


112 


C 


264 


04 54 38.54 


+03 00 51.5 


17.43 


61074 


104 


C 


265 


04 54 39.39 


+03 03 44.7 


20.44 


60825 


189 


C 


266 


04 54 40.49 


+02 52 31.8 


19.67 


66941 


144 


C 


267 


04 54 40.76 


+03 02 35.9 


19.51 


60006 


117 


C 


268 


04 54 41.11 


+02 59 49.5 


21.97 


113621 


450 


C 


269 


04 54 41.67 


+02 59 17.8 


20.89 


106558 


108 


C 


270 


04 54 41.90 


+02 58 55.2 


19.60 


69483 


270 


C 


271 


04 54 42.25 


+03 02 02.1 


20.33 


99066 


126 


C 


272 


04 54 42.63 


+03 01 58.3 


19.24 


59602 


144 


C 


273 


04 54 43.01 


+02 57 42.3 


20.60 


70781 


72 


C 


274 


04 54 44.96 


+03 01 06.2 


20.78 


113447 


225 


C 


275 


04 54 45.34 


+02 53 40.4 


19.59 


55623 


104 


C 


276 


04 54 46.28 


+02 52 36.1 


20.56 


77385 


99 


C 


277 


04 54 46.85 


+02 53 21.9 


20.84 


55396 


99 


C 


278 


04 54 46.85 


+03 02 57.0 


21.49 


98017 


99 


C 


279 


04 54 47.88 


+03 03 33.3 


18.62 


71557 


112 


c 


280 


04 54 48.86 


+02 52 34.5 


20.66 


62366 


180 


c 



Table 1. Continued. 
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281 


04 54 49.21, +03 03 23.4 


20.20 


99588 


180 


C 


282 


04 54 51.40, +02 57 52.4 


19.39 


45940 


99 


C 


283 


04 54 52.12, +02 54 27.4 


20.81 


69405 


117 


c 


284 


04 54 52.32, +02 56 42.7 


20.45 


110251 


225 


c 


285 


04 54 53.87, +02 52 45.2 


20.03 


60537 


126 


c 


286 


04 54 55.00, +02 54 00.0 


19.21 


62489 


148 


c 


287 


04 54 55.47, +02 58 31.7 


22.45 


138135 


144 


c 


288 


04 54 56.66, +03 00 53.4 


21.63 


94455 


112 


c 


289 


04 54 56.87, +02 54 02.9 


18.69 


69321 


99 


c 


290 


04 54 57.20, +03 03 26.6 


20.64 


98625 


104 


c 


291 


04 54 58.37, +03 01 54.0 


20.73 


135512 


225 


c 


292 


04 54 59.89, +02 53 31.1 


19.46 


68559 


135 


c 


293 


04 55 01.66, +02 57 55.7 




138036 


225 


c 



center of the whole cluster, the optical (biweight) center of the 
HVG is closer to the X-ray cluster center (see Fig. |7]i. 
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Fig. 3. Spatial distribution on the sky of the 293 galaxies hav- 
ing redshifts in the cluster field. Circles indicate galaxies having 
new redshifts acquired with the TNG. Dots and triangles indi- 
cate galaxies having redshift data from CNOC and Proust et al. 
(120001) catalogs, respectively. The X-ray peak is taken as the 
cluster center. 
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Fig, 4. Redshift galaxy distribution. The solid line histogram 
refers to the (223) galaxies assigned to the cluster according to 
the DEDICA reconstruction method. 



We also use the results of the gap analysis to determine the 
first guess when using the Kaye's mixture model (KMM) test 
to find a possible group partition of the velocity distribution (as 
implemented by Ashman et al. |1994| l. The KMM algorithm fits 
a user-specified number of Gaussian distributions to a dataset 
and assesses the improvement of that fit over a single Gaussian. 
In addition, it provides the maximum-likelihood estimate of the 
unknown n-mode Gaussians and an assignment of objects into 
groups. We find a two-groups partition which is a significantly 
better descriptor of the velocity distribution with respect to a sin- 
gle Gaussian at the 95% cl.. The cluster partition is similar to 
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Fig. 5. Lower panel: projected clustercentric distance vs. rest- 
frame velocity for the 223 galaxies in the main peak (Fig. H]). 
Crosses show galaxies detected as interlopers by our "shifting 
gapper" procedure. Upper panel: rest-frame velocity histogram 
for the 223 galaxies in the main peak; the solid line refers to 
the 167 cluster members only. Large and small arrows indicate 
the positions of weighted gaps in the velocity distribution of the 
whole sample and of the main system (MS). Labels bl and b2 
indicate backl and back2 "background" peaks of galaxies. 

Table 2. Results of the weighted gap analysis for the whole sam- 
ple and for the MS subsystem. 



Fig. 6. Integral profiles of mean velocity {upper panel) and LOS 
velocity dispersion {lower panel). The mean and dispersion at a 
given (projected) radius from the cluster center is estimated by 
considering all galaxies within that radius (the first point is ob- 
tained on the basis of the five galaxies close to the cluster center). 
The error bands at the 68% c.l. are shown. In the lower panel, 
the horizontal line represents the X-ray temperature with the re- 
spective 90 per cent errors (Govoni et al. 2004) transformed in 
CTv assuming the density-energy equipartition between gas and 
galaxies, i.e. ySspec - 1 (see text). 

Table 3. Global properties of the whole sample, the MS and the 
HVG. 



Sample 


Ngaljj pre>Ngais aft 


Vpre, Vaft 


Size 


Prob. 








km s"' 








Whole Sample 


145,22 


61547,61964 


3.81 


5.QE- 


4 


MS 


30, 30 


59059,59163 


2.52 


lAE- 


2 


MS 


30,33 


59722, 59783 


2.35 


3.0E- 


2 


MS 


33,52 


60440,61548 


2.34 


3.0E- 


2 





Sample 


N„ 


< V > 
kin s"' 


CTv 

km s"' 


Rvii 
h-^'Mpc 


Mass(< Rvir) 


Whole system 

MS 

HVG 


167 
145 
22 


60209 + 82 
59978 + 67 
62419 ± 74 


1066:« 

338+^f 


2.34 
1.79 
0.74 


17 ±2 
8 + 2 

6+" ** 

-0.3 



that indicated by the above weighted gap analysis detecting two 
groups with 146 and 21 galaxies. 

4.3.2. Dressler-Shectman statistics 

We also analyze substructure combining galaxy velocity and 
position information. We compute the A-statistics devised by 
Dressier & Shectman (1988 , hereafter DS). We find a significant 
indication of DS substructure (at the 97% c.l. using 1000 Monte 
Carlo simulations; see e.g. Boschin et al. 2004). Figure [8] shows 
the distribution on the sky of all galaxies, each marked by a cir- 
cle: the larger the circle, the larger the deviation <5, of the local 
kinematical parameters from the global cluster parameters, i.e. 
the higher the evidence for substructure. 

To better point out galaxies belonging to substructures, we 
resort to the technique developed by Biviano et al. ( 120021 see 



also Boschin et al. 120061 Girardi et al. I2006I I. who used the indi- 
vidual 5,-values of the DS method. The critical point is to deter- 
mine the value of that optimally indicates galaxies belonging 
to substructure. To this aim we consider the i5,-values of all 1000 
Monte Carlo simulations used above. The resulting distribution 
of Si is compared to the observed one finding a difference at the 
99% c.l. according to the IDKS-test. The "simulated" distribu- 
tion is normalized to produce the observed number of galaxies 
and compared to the observed distribution in Fig. |9] the latter 
shows a tail at large values. Selecting galaxies with 6i < 2.5 the 
IDKS-test gives only a marginal difference between real and 
simulated galaxies (at the 93% c.l.) suggesting that galaxies with 
6i > 2.5 presumably are in substructures. These galaxies are in- 
dicated with heavy circles in Fig. [8] showing four subclumps (at 
northern, eastern, southern and distant western cluster regions). 
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Fig. 7. Spatial distribution on the sky of the 167 galaxies of the 
whole cluster showing the two groups recovered by the weighted 
gap analysis. Dots and small circles indicate the main system 
(MS) and the high velocity group (HVG) galaxies. Large cross, 
rotated square and circle indicate the optical centers of the whole 
cluster, the MS and the HVG, respectively. The X-ray peak is 
taken as the cluster center 
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Fig. 8. Spatial distribution of the 167 cluster members, each 
marked by a circle: the larger the circle, the larger is the devi- 
ation 5i of the local parameters from the global cluster param- 
eters, i.e. there is more evidence for substructure (according to 
the Dressier & Shectman test, see text). Heavy circles indicate 
those with (5, > 2.5. 



4.3.3. Analysis of velocity dispersion profiles 

Finally we analyze the kinematical properties of galaxy pop- 
ulations located in different spatial regions of the cluster We 
compute the profiles of mean velocity and velocity dispersion 
of galaxy systems surrounding the lensing mass peaks listed 
by M07 (see Fig. [TO] i. This allows an independent analysis of 
the possible individual galaxy clumps. A quasi flat profile is 
expected in the case of a relaxed system with isotropic orbits 
for galaxies (e.g., Girardi et al. 1998). Although an increas- 
ing/decreasing profile might be due to particular orbits of galax- 
ies in a relaxed system (e.g., Girardi et al. 1998; Biviano & 
Katgert 2004), here this is likely connected with the presence 
of substructure. As for an increasing profile, this might be sim- 
ply induced by the contamination of the galaxies of a close, sec- 




Fig. 9. The distribution of 6/ deviations of the Dressler- 
Shectman analysis for the 167 member galaxies. The solid line 
represents the observations, the dashed line the distribution for 
the galaxies of simulated clusters, normahzed to the observed 
number. 



ondary clump having a different mean velocity (e.g., Girardi et 
al. 119961 Girardi et al. .2006) . This hypothesis can be investi- 
gated by looking at the behavior of the mean velocity profile. 
In fact, if the o-y profile increases due to the contamination of a 
close clump, for the same reason and at about the same radius, 
the < v > profile should increase/decrease. As for a decreasing 
profile, this might be likely due to the projection effect of a few 
clumps centered around the center of the system and having dif- 
ferent mean velocities, i.e. somewhat aligned with the LOS or, 
alternatively, of a large scale structure (LLS) elongated along the 
LOS (e.g. a LLS filament). 

The inspection of Fig.[TO]shows that the o-y profiles of peaks 
No. 1, 2 and 5 sharply increase with the distance from the peak 
position. Simultaneously, the < v > profiles decline (peak No. 
1) or increase (peaks No. 2 and 5). For each of these clumps 
we attempt to detect the region likely not contaminated by other 
clumps - and thus reliable for kinematical analysis - as the re- 
gion before the sharp increasing of the o-y profile (see the arrows 
in Fig.fTOlfor peaks No. 1, 2, and 5). No conclusion can be driven 
for peaks No. 3 and 4 where the cry profile is decreasing. 

4.4. Substructure of the main system 

Here we present the results of our substructure analyses applied 
to the 145 galaxies of the main system (MS), i.e. rejecting the 
galaxies of the high velocity group (HVG) which might mask 
the real cluster structure. 

Table |4] and Figure [TT| summarize kinematical and spatial 
properties of the subclumps we detect in the MS. In particular, 
Table|4]lists for each clump the number of galaxies (Col. 2); the 
mean velocity and its jacknife error (Col. 3); the velocity dis- 
persion and its bootstrap error (Col. 4); the luminous galaxies 
contained within the analyzed clump (Col. 5); the name of the 
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Fig. 10. Kinematical profiles of the galaxy clump surrounding 
the lensing mass peaks No. 1-5 listed by M07. For each peak in- 
tegral mean velocity (in units of 10^ km s"') and LOS velocity 
dispersion (in units of km s ') profiles are shown in upper and 
lower panels, respectively. The technique is the same adopted 
in Fig. |6l but we omit errors for the sake of clarity. The solid 
and dashed lines join the values obtained using all cluster galax- 
ies and only galaxies belonging to the MS, respectively. The 
arrows indicate the regions likely not contaminated from other 
subclumps (see Sect. 14.3. 3l and Sect. 14.4.31 ). 



corresponding structure discussed in Sect. |5] (Col. 6). The fol- 
lowing subsections show the results recovered for each of the 
three methods of analysis. 
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Fig. 11. Summary of spatial distribution of cluster substruc- 
ture. Labels "MS" and "HVG" indicate the optical centers of 
the MS and of the HVG, respectively. Other symbols refer to 
galaxy subclumps found in the MS using several, different ap- 
proaches. Triangles, squares, circles and rotated squares indi- 
cate 2D DEDICA peaks (large symbols) and luminous galaxies 
(small symbols) of the VI, V2, V3 and V4 clumps. Labels "N", 
"S", "E" and "W" indicate the (biweight) centers of the four DS 
clumps. Labels "1", "2", "3", "4" and "5" indicate the corre- 
sponding lensing mass peaks listed by M07: around peaks No. 
1, 2, 4 and 5 we detect the PI, P2, P4 and P5 clumps. On the 
base of the location of these subclumps and their velocities in 
Table |4] we discuss the presence of five main structures {N&l, 
N&2, S'W, & and 'W centered around the labels "1", "2", "4", 
"5" and "W", respectively) and of two minor structures (N and 
C roughly located around the label "N" and the large rotated 
square, respectively), see Sect. |5] The arrow indicates the head 
tail radiogalaxy ID 184 with the direction of its tail. The plot is 
centered on the X-ray cluster center. 



4.4.1. Velocity distribution 

The velocity distribution of the MS is negatively skewed (at the 
c.l. of 90-95%, skewness=-0.330) and Ught-tailed (at the c.l. of 
90 - 95%, kurtosis=2.347). The W-test (Shapiro & Wilk [T965] l 
rejects the null hypothesis of a Gaussian parent distribution at 
the > 99.9% c.l.. 

We detect three marginally significant gaps which divide the 
MS in four groups of 30, 30, 33 and 52 galaxies (see Fig.|5]and 
Table |2]), hereafter defined as VI, V2, V3 and V4 from low to 
high velocities. When compared two by two through the 2DKS- 
test, these groups differ in spatial distribution: the VI group dif- 
fers both from the V2 and V3 groups (at the 98% and 94% c.l., 
respectively); the V4 group differs from the V3 group (at the 
94% c.l.). For each of the these groups Fig.[T2lshows the spatial 
distribution of galaxies and the corresponding peak according to 
the 2D DEDICA procedure. For the V2 group, we find and plot 
two peaks of comparable significance. Figure [TT] shows the po- 
sition of these peaks in relation to the lensing mass peaks listed 
by M07. Each of these groups contains one o more luminous 
galaxies and, with the exception of V4, each group hosts one 
correspondent luminous galaxy close to the respective peak in 
the galaxy distribution (see Fig.fTTTi. Properties of groups recov- 
ered by kinematical analysis are listed in Table H] For the VI, 
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Fig. 12. Spatial distribution on the sky of the 145 galaxies be- 
longing to the MS showing the four groups recovered by the 
weighted gap analysis. Small triangles, squares, circles and dots 
indicate galaxies of the VI, V2, V3 and V4 groups, respectively. 
Large triangle, (2) squares, circle and rotated square indicate 
DEDICA peaks of the VI, V2, V3 and V4 galaxy distributions. 



M. Girardi et al.: Cluster Abell 520: a 
Table 4. Results of the substructure analysis of the MS. 



Gal. Clump N„ < v > trj Lum.gals Structures 

km s"' km s ' ID 

VI'' 30 58568 ±56 300!^^ 205 £ 

¥2" 30 59455 ± 18 95!* 44,170 N&2 + ^V 

V3'' 33 601 13 ±19 106+1^ 204, 160, 137 ATfil 

V4'' 52 60920 ±33 238+' 106,264 S'W + C 

DS-N 4 58427 ±272 415!*|* - 

DS-S 4 61165 ±519 19\*f}^ 106 S'W 

DS-E 7 58686 ±292 682+^56 205 £ 

DS-W 4 59635 ±249 380+*3, - ^V 

PI 7 60448 ±349 811^^™ 204 N8\ 

P2 9 59596 ±276 749^'°* 170,160 A^£2 

P4 6 61447 ±283 579:=| 106 S'W 

P5 6 58634 ±325 668+^™ 205 £ 



We use the biweight and the gapper estimators by Beers et al. 
(1990) for samples with Ng > 15 and with Ng < 15 galaxies, re- 
spectively (see also Girardi et al. |l993i l. 

The estimate of cry should be considered a lower limit in these sam- 
ples (see text). 



V2, V3 and V4 groups, the membership assignment might lead 
to an artificial truncation of the tails of the distributions; thus the 
values of velocity dispersion should be considered lower limits 
(e.g., Bird fT994l) . 

Using the results of the gap analysis to determine the first 
guess of the KMM algorithm we find that a four-groups partition 
is a significant better descriptor of the velocity distribution with 
respect to a single Gaussian at the 99% cl.. In particular, the clus- 
ter partition is similar to that indicated by the above weighted- 
gap analysis separating the MS in groups of 33, 27, 35 and 50 
galaxies. 

4.4.2. Dressler-Shectman statistics 

The DS test on the MS gives a very marginal indication of sub- 
structure (at the 91% c.l.). However, as for the location of sub- 
structures, the DS plot of the MS is similar to that recovered for 
the whole sample (cf. Fig. [8] and Fig. [T3] l. The main difference 
is the disappearance of the northern substructure very close to 
the X-ray cluster center which was likely due to the galaxies of 
the HVG. The other three subclumps are still present (at east- 
ern, southern and distant western cluster regions). Moreover, the 
analysis of the MS shows a northern subclump 5' from the X- 
ray cluster center. In order to better investigate the properties of 
DS substructure we select the galaxies with the highest 5, in such 
a way that at least four galaxies can be assigned in an unique way 
to each of the four DS subclumps. This leads to 19 galaxies with 
(5, > 1.95. The kinematical properties of the four DS subclumps 
are listed in Table H] Figure [TT] shows the position of their (bi- 
weight) centers. 

4.4.3. Analysis of velocity dispersion profiles 

We also reanalyze the kinematical properties of galaxies sur- 
rounding the tensing mass peaks listed by M07 (see Fig. [10] 
- dashed lines). We notice that for all the peaks where we 
have defined a region likely uncontaminated by other clumps in 
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Fig. 13. The result of Dressler-Shectman analysis as in Fig. [8] 
but for the 145 members of the MS. Here heavy circles indicate 
galaxies with 5, > 1.95 used to define the DS-N, DS-S, DS-E 
and DS-W clumps (see text). 

Sect. |4.3.3l (peaks No. I, 2 and 5), we confirm that these regions 
are devoid of HVG galaxies. Moreover, now we also find a sharp 
increase of the cr^ profile for peak No. 4 and we define a likely 
non contaminated region for the subsystem corresponding to this 
peak, too. The groups formed by galaxies within the uncontami- 
nated fiducial regions are referred as PI, P2, P4 and P5 and their 
kinematical properties are shown in Table|4] 

In the case of peak No. 3 we still find a decreasing cr^ profile; 
thus we find no evidence for an individual, dynamically impor- 
tant structure around this peak. 

5. Cluster structure and dynamics 

The value we find for the global velocity dispersion of the cluster 
members, cry — I066+gj, is in agreement with previous analyses 
(Carlberg et al. fT996l Borgani et al. 1999; Proust et al. 20001 
Mezzetti & Girardi 2001). This value of the velocity dispersion 
is also comparable to the average X-ray temperature assuming 
the density-energy equipartition between gas and galaxies, i.e. 
Apec = 1 B see Fig.|6]- lower panel. 

We analyze the cluster structure using the velocity distribu- 
tion analysis (weighted gap technique and KMM method), the 
Dressler-Shectman statistics and the analysis of the velocity dis- 
persion profiles. The structure of A520 is definitely very com- 
plex and thus likely far from the dynamical equilibrium. The 
agreement between cr^ and Tx might be due to an enhancement 
of both quantities and, as already pointed out by M07, gross 
properties and scaling relations are not always useful indicators 
of the dynamical state of clusters. 

Hereafter we analyze and discuss the structure of A520 start- 
ing from the simplest hypothesis and then adding some degrees 
of complexity. 

5.1. Mass estimate 

Making the usual assumptions (cluster sphericity, dynamical 
equilibrium, that the galaxy distribution traces the mass distri- 
bution), one can compute virial global quantities. Following the 

3 jSspec = / (kT / ixnip) with /i = 0.58 the mean molecular weight and 
lUp the proton mass. 
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prescriptions of Girardi & Mezzetti (1200 11 1, we assume for the 
radius of the quasi-virialized region R^h = 0.17 x crylH(z) - 
2.34 /i^d Mpc - see their eq. 1 after introducing the scaling with 
Hiz) (see also eq. 8 of Carlberg et al. [19971 for R2im). Thus the 
cluster is sampled out - although in a non-homogeneous way - 
to 7?200- We compute the virial mass (Limber & Mathews 1 19601 
see also, e.g., Girardi et al. |l998] i: 

M = 37r/2 ■ o-IRvm /G - SPT, ( 1 ) 

where SPT is the surface pressure term correction (The & White 
1 19861 1. and ^pv is a projected radius (equal to two times the pro- 
jected harmonic radius). 

The estimate of cry is robust when computed within a large 
cluster region (see Fig.|6l). The value of 7?pv depends on the size 
of the sampled region and possibly on the quality of the spa- 
tial sampling (e.g., whether the cluster is uniformly sampled or 
not). Considering the 155 galaxies within Ry„ we obtain /?pv - 
(1.72 ± 0.09) h^iQ Mpc, where the error is obtained via the jack- 
nife procedure. The value of SPT strongly depends on the radial 
component of the velocity dispersion at the radius of the sampled 
region and could be obtained by analyzing the (differential) ve- 
locity dispersion profile, although this procedure would require 
several hundred galaxies. We decide to assume a 20% SPT cor- 
rection as obtained in the literature by combining data on many 
clusters sampled out to about /?vir (Carlberg et al. il997i: Girardi 
et al. |T998] l. We compute M(< R^,, = 2.34 /i^ijMpc) = (1.7±0.2) 

xlO^^ Mq. 

Since the cluster center of A520 is not well defined and the 
spatial sampling is not complete and homogeneous within 7?vir, 
one could use an alternative estimate of /?pv on the basis of the 
knowledge of the galaxy distribution. Following Girardi et al. 
( 1199 8 ; see also Girardi & Mezzetti l2001l) we assume a King-like 
distribution with parameters typical of nearby/medium-redshift 
clusters: a core radius R^ - 1/20 x /?vir and a slope-parameter 
/?fit — 0.8, i.e. the volume galaxy density at large radii goes as 
= r~2-4. We obtain R^m = 1.74 li'l Mpc, where a 25% 
error is expected (Girardi et al. 11998 ). The mass recovered by 
this method is then M(< R^i, = 2.34 /j^o'Mpc) = (1.7 + 0.5) 
xlO'^ h~Q Mq in excellent agreement with the above direct esti- 
mate. 

Our analysis of the cluster velocity distribution detects the 
presence of a high velocity group (HGV) with a relative rest- 
frame LOS velocity of Vrf ~ 2000 km s"' with respect to the 
main system (MS), see § 14.3.11 Therefore we might think that 
the cluster is better described by the combination of the MS 
and the HVG, considered as two separated entities. Assuming 
the dynamical equilibrium for both the MS and the HVG we 
can compute independent virial radii and masses M(< 7?vii - 
1.79 /i^o'Mpc) = (8 + 2) xlO'"* h'^ Mq and M{< R^i, = 
0.74 h-^Mpc) = (0.6+°-^) xlO'-* /j^,,' Mq, respectively. 

To compare the mass estimates derived for the two above 
cluster models (relaxed cluster and the MS-hHVG system) we 
consider the mass values within 1 /i^g Mpc . To rescale our 
mass estimates we assume that the system is described by a 
King-like mass distribution (see above) or, alternatively, a NFW 
profile where the mass-dependent concentration parameter c is 
taken from Navarro et al. ( 1997) and rescaled by the factor 1 + z 
(Bullock et al. l200n Dolag et al. l2004] l. i.e. c = 4.16, 4.96 and 
6.18 for the whole cluster, the MS and the HVG, respectively. 
The relaxed cluster model leads to a mass of M(< 1 /i^gMpc) - 
(7.2 - 9.6)xl0''* h-^ Mq while the addition of the MS and the 
HVG masses leads to a mass of M(< 1 /i^^jMpc) = (4.0 - 6.8) 



xlO'"^ /i^Q Mq [where the mass range includes a Icr error on the 
original M{< /?vii-) estimate]. 

Using the above rescaling we can also compare our results 
with the estimates recovered from X-ray and gravitational lens- 
ing analyses. Lewis et al. (,1999) used the ROSAT X-ray sur- 
face brightness and ASCA temperature to estimate a mass of 
M(< 1.764 h-^Mpc) = (11.3 + 1.1) xlO'"* /i^,] Mq (the 
ROSAT-PSPC estimate is converted in our cosmology). X-ray 
mass is intermediate between our estimates since we obtain 
M(< 1.764 /!7(jMpc) = (12 - 15) xlO''* h'^ Mq and (6.1 - 10) 
xlO''* /z^g Mq for the two cluster models. As for gravitational 
lensing, Dahle et al. (I2002I I obtained the projected mass Mproj(< 
1.111 /i^o'Mpc) = n.7+^-^ xlO'"* Mq (see their Fig. 50 with 
conversion in our cosmology). Our projected mass estimates for 
the two cluster models are Mproj(< 1.111 /i^ijMpc) = (12 - 18) 
xlO'"* Mq and (5.4 - 1 1) xlO'"* h'l^ Mq, where to make the 
projection we have considered that the cluster mass distribution 
is truncated at one or at two virial radii. 

5.2. Main system and high velocity group: relative dynamics 

Continuing with the assumption of a cluster formed by the MS 
and the HVG, we investigate their relative dynamics. We use 
different analytic approaches which are based on an energy in- 
tegral formalism in the framework of locally flat spacetime and 
Newtonian gravity (e.g.. Beers et al. 119821 ). The values of the rel- 
evant observable quantities for the two-clumps system are: the 
relative LOS velocity in the rest frame, Vrf = 2033 km s"' (as 
recovered from the MS and the HVG); the projected linear dis- 
tance between the two clumps, D = 0.21 /i^,] Mpc (as recov- 
ered from optical centers of the MS and the HVG); the mass of 
the system obtained by adding the masses of the two subclus- 
ters each within its virial radius, logMjys - 14.9154;^g [gj2 (see 
Tableig. 

First, we consider the Newtonian criterion for gravita- 
tional binding stated in terms of the observables as V^D < 
IGMsysSin^acosa, where a is the projection angle between the 
plane of the sky and the line connecting the centers of the two 
clumps. The thin curve in Fig. [14] separates the bound and un- 
bound regions according to the Newtonian criterion (above and 
below the curve, respectively). Considering the value of Msys, 
the MS+HVG system is bound between 21° and 83°; the cor- 
responding probability, computed considering the solid angles 

(i.e., J^j cos a da), is 63%. 

Then, we apply the analytical two-body model introduced 
by Beers et al. (|1982| l and Thompson (119821 see also Lubin et 
al. 1998 for a recent application). This model assumes radial 
orbits for the clumps with no shear or net rotation of the sys- 
tem. Furthermore, the clumps are assumed to start their evolu- 
tion at time fo = with separation do = 0, and are moving 
apart or coming together for the first time in their history; i.e. 
we are assuming that we are seeing the HVG prior to merging 
with the MS (at the time t= 1 1 .022 Gyr at the cluster redshift, see 
Wright [2006l l. The bimodal model solution gives the total system 
mass Msys as a function of a (e.g., Gregory & Thompson |1984| l. 
Figure [T4| compares the bimodal-model solutions with the ob- 
served mass of the system. The present bound outgoing solutions 
(i.e. expanding), BO, are clearly inconsistent with the observed 
mass. The possible solutions span these cases: the bound and 
present incoming solution (i.e. collapsing), Bla and Bib, and the 
unbound-outgoing solution, UO. For the incoming case there 
are two solutions because of the ambiguity in the projection an- 
gle a. We compute the probabilities associated to each solution 
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Fig. 14. System mass vs. projection angle for bound and un- 
bound solutions (solid and dotted curves, respectively) of the 
two-body model applied to the MS and the HVG subsystems. 
Labels Bla and Bib indicate the bound and incoming, i.e. col- 
lapsing solutions (the main part of the solid curve). Label BO 
indicates the bound outgoing, i.e. expanding solutions (the part 
of the solid curve which is roughly a vertical line). Label UO in- 
dicates the unbound outgoing solutions (the dotted curve which 
is roughly a vertical line). The horizontal lines give the observa- 
tional values of the mass system and its uncertainties. The bound 
and unbound regions according to the Newtonian criterion are 
indicated, too (above and below the thin curve, respectively). 



assuming that the region of Msys values between the uncertain- 
ties are equally probable for individual solutions: f Bia ~ 88%, 
f Bib ~ 12%, Puo ~ 7 X 10"^%. Thus it is very likely that the 
HVG lies in front of the cluster just infalling onto it. 

Notice, however, that the centers of the HVG and the MS 
are not well determined. The HVG has too small a number of 
galaxies for a precise center determination. As for the MS center, 
we might adopt the X-ray cluster center instead of the optical 
cluster center These uncertainties do not change the bulk of our 
results. For instance, in the case where we assume that the X-ray 
center is the MS center, i.e. a smaller projected linear distance 
between the two clumps {D = 0.1 /j^q Mpc), the effect is to 
increase the boundary probability (at 75% c.l. for the Newtonian 
model) and to yield more extreme values for the bound solutions 
(a ~ 15 and 85 degrees). Also, possible underestimates of the 
masses (e.g., if the MS and the HVG actually extend outside 
of the virial radii we estimate for them) would lead to binding 
probabilities larger than those computed above, as well as more 
extreme values for a. Thus the analysis here displayed should be 
considered a lower limit for our conclusions in § 15.51 where we 
propose the existence of a cluster accretion along the LOS (i.e. 
a close to 90 degrees). 



5.3. NE-SW merger 

Although the presence of the HVG system is maybe the most 
important for the optical virial mass computation, we find that 
A520 shows a much more complex structure. In fact, we detect 
several subsystems in the MS along the NE-SW and E-W direc- 
tions (see Figs. [Tl and fTTTi. 

Clumps around lensing mass peaks No. 1, 2 and 4 of M07 
define the direction of the likely merger along the NE-SW di- 
rection with the SW structure having crossed the NE structure 
(Markevitch et al. l2005] l. 

Around peaks No. 1 and 2 we detect two structures (here- 
after N&l and N&2) using the analysis of the velocity distri- 
bution (V3 and V2) and velocity dispersion profiles (PI and 
P2). These structures have a relative rest-frame LOS velocity 
of ~ 700 km s '. Each structure hosts, close to the center, a lu- 
minous galaxy having typical velocity of the structure (IDs 204 
and 170 in N&\ and N&2, respectively). 

Structures N&\ and N&2 are so close in space and in ve- 
locity that the "uncontaminated" regions detected by cr^ profiles 
are slightly superimposed (see Fig.[Tli. Moreover, at < 0.5 'from 
ID 170 (the N&2 central, luminous galaxy) we find the luminous 
galaxy ID 160 which belongs to V3 group, i.e. to the N&\ struc- 
ture. Such close couples of galaxies having different velocities 
are often observed in clusters (Boschin et al. |2006I Barrena et 
al. 2007a) and are the likely tracers of a previous cluster merger. 
Indeed cluster merger is thought to be the cause of the forma- 
tion of dumbbell galaxies (e.g.. Beers et al. 119921 Flores et al. 
I2000I I. Therefore N&1+N&2 is likely to form a single, although 
not yet well relaxed, structure and represent the real, original 
main cluster. In fact, the combined velocity of N&l and N&2 
is ~ 60000 km s"', similar to the mean velocity of the MS, and 
the N&2 position is close to the optical center of the MS. The 
dynamical importance of the N&\+N&2 structure explains why 
the SW structure has been reported to have suffered significant 
damage in the merging, as shown by the pieces of the cluster 
core detected in X-ray (Markevitch et al. 12005b . 

Around peak No. 4 we detect a structure (hereafter S'W) 
both using DS analysis (DS-S) and studying the velocity disper- 
sion profile (P4). It is characterized by a high velocity v ~ 6 1 300 
km s"', i.e. v,f ~ +1100 km s"' from the N&1+N&2 complex. 
S'W is not individually detected in the velocity distribution. 
However, since it hosts close to its center one of the two lumi- 
nous galaxies of the V4 group (ID 106 with v = 61277 km s"'), 
S'W is likely a part of V4 (see also discussion in Sect. l5.5l l. 

A possible structure related to the NE-SW merger is DS-N 
(hereafter AO having a small velocity v ~ 58400 km s"'. It is 
only detected through the Dressler-Shectman analysis and does 
not contain a luminous galaxy; therefore we neglect it in the fol- 
lowing discussion. Notice however that it roughly corresponds to 
the N peak in the lensing mass map of Okabe & Umetsu (20081 1. 

To attempt a more detailed analysis of the NE-SW merger 
we apply the bimodal model considering the interaction between 
N&1+N&2 (likely corresponding to the main part of the MS, 
see above) and S'W assuming we are looking at them after their 
core crossing as suggested from X-ray data. As parameters of 
the model we use a relative LOS velocity V^ = 1 100 km s ' and 
a relative projected distance D - 0.8 /i^^ Mpc . We assume the 
MS mass as the mass for the whole system. In order to apply the 
two-body model we assume that the time to -0 with separation 
(io = is the time of their core crossing and that we are seeing 
the cluster a few 10^ years after the merging. In fact, a few 10^ 
years is the time scale in which the relativistic electrons lose en- 
ergy, i.e. the lifetime of the radio halo (e.g., Giovannini & Feretti 
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Fig. 15. System mass vs. projection angle for bound outgoing 
and unbound solutions (solid and dotted curves, respectively) of 
the two-body model applied to the N&1+N&2 and S'W sub- 
systems. Labels Blt2, BOt2 and UOt2 indicate the curves cor- 
responding to bound incoming, bound outgoing and unbound 
outgoing solutions assuming that the system is observed at the 
time t -0.2 Gys after the core crossing. Labels Blt3, B0t3 and 
UOt3 have the same meaning but assuming that the system is 
observed at t =0.3 Gys after the core crossing. The horizontal 
lines give the observational values of the mass system and its 
uncertainties. The bound and unbound regions according to the 
Newtonian criterion are indicated, too (above and below the thin 
curve, respectively). 



12002 ). Figure [T5] shows the results for a time of f =0.2 and 0.3 
Gyrs after the core crossing. The likely bound solution is then an 
outgoing one, i.e. the SW structure is now behind the NE struc- 
ture going away from it. In particular, an angle of or ~ 30 leads 
to an outgoing velocity of ~ 2200 km s"'. 

The Mach number of the shock is Al = Vj/cs, where Vj is 
the velocity of the shock and Cs in the sound speed in the pre- 
shock gas (see e.g., Sarazin |2002| for a review). In the station- 
ary regime we can assume that Vs is the merger velocity 2200 
km s"'. Assuming the equipartition of energy density between 
gas and galaxies and using the o-y of the MS we obtain ~ 
812 km s ' from the thermal velocity. Therefore we estimate 
M ~ 2.7, which is in reasonable agreement with M - 2.1;^q3 
and M = 2.2^y j recovered from X-ray data (Markevitch et al. 
l2005T l. 

Notice that f = 0.2 - 0.3 Gyrs and o- ~ 30 degrees are 
larger, but comparable, to the merger parameters recovered by 
Markevitch et al. (2002) for a similar scenario in the cluster 
1E0657-56 which hosts a radio halo and shows a bow shock, 
too (0.1-0.2 Gyrs and 10-15 degrees). Indeed, this agreement is 
not casual since as noticed by Markevitch et al. (I2005I I to ob- 
serve a shock front "one has to catch a merger at a very specific 
stage when the shock has not yet moved to the outer, low surface 
brightness regions and at a sufficiently small angle from the sky 
plane, so that projection does not hide the density edge". 



5.4. Accretion along the E-W direction 

M07 suggested a possible secondary E-W merger related to 
lensing mass peaks No. 3 and 5. However, these authors doubted 
the dynamical importance of the structure around the peak No. 5 
since this mass peak is poorly significant and the mass-to-light 
ratio is quite low (but see Okabe & Umetsu l2008i where the cor- 
responding C3 peak is quite significant). 

We find a strong dynamical evidence of a structure around 
peak No. 5 (hereafter £) at v ~ 58600 km s"' (i.e. Vrf ~ -1 150 
km s"' from the N&1+N&2 complex) using the velocity distri- 
bution analysis (VI), the velocity dispersion profile (P5) and DS 
analysis (DS-E). Moreover, & hosts its luminous galaxy ID 205 
close to its center Comparing & with S'W we notice that & lies 
at a similar velocity distance and at much smaller spatial dis- 
tance with respect to the N&1+N&2 complex. However, since 
there is no sign of a strong, present interaction from X-ray data, 
we suspect that £ might be a high-speed remnant of a previous 
merger. 

We also detect a structure (hereafter "W) located in the west- 
ern external region not sampled by previous gravitational lensing 
analyses. This structure is found using both the velocity distri- 
bution analysis (part of V2) and the DS analysis (DS-W) and is 

characterized by a velocity v ~ 59600 km s"' (i.e. Vrf 300 

km s"' from the N&\+N&2 complex). It hosts a luminous 
galaxy close to its center, too (ID 44). Since the velocity is 
similar to that of the N&\ +N&2 complex and the projected spa- 
tial distance is about 2 Mpc, might be a distant subclump 
well far from the merging, infall phase. 

In conclusion we strongly reinforce the possibility of an ac- 
cretion onto A520 along the E-W direction. 



5.5. Nature otthe "massive dark core" 

Finally we discuss the region around lensing mass peak No. 3 for 
which M07 found a very large mass-to-light ratio claiming for 
the presence of a massive dark core (but see Okabe & Umetsu 
120081 where the corresponding CI peak is not particularly pro- 
nounced). 

This peak is the only M07 peak for which our analysis of the 
velocity dispersion profile does not support the presence of an 
individual structure. The peak of the V4 galaxy distribution is 
close to peak No. 3, but the V4 group is the only one which does 
not host any luminous galaxy close to its center In fact, the two 
V4 luminous galaxies are located well far in the northern and 
southern cluster regions (see Fig. fTTT l. Moreover, several galax- 
ies of V4 are likely to be associated with the V4 southern lumi- 
nous galaxy, i.e. with the S'W structure we discuss in Sect. 15.31 
Therefore, only part of galaxies we assign to V4 are likely really 
connected with the region around peak No. 3 (hereafter we name 
this minor structure C). 

In conclusion, the existence of an individual, very important 
structure associated to peak No. 3 is not supported by our kine- 
matical analysis. Rather, we find evidence of two groups cen- 
tered in that region: the minor C clump (at v~ 60900 km s ') 
and the HVG (at v~ 62400 km s"') suggesting the accretion 
onto the cluster along the LOS. In agreement with this idea, in 
the MS+HVG bimodal model of Sect. 15.21 we prefer the bound 
solution with a > 80 degrees, with the HVG almost LOS aligned 
with the MS and infalling onto it. A scenario of a few groups at 
different velocities agrees with the high velocity dispersion of 
galaxies we measure around this peak, see also the high velocity 
ridge in the merging cluster Abell 521 (Ferrari et al. i2003i ). 
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Fig. 16. Spatial distribution on the sky of the 293 galaxies hav- 
ing redshifts in the cluster field. Dots and large symbols indi- 
cate member and non member galaxies, respectively. In particu- 
lar, squares and circles indicate galaxies belonging to backl and 
back2 groups. 



Working in this scenario we also look for other possible 
groups aligned with the LOS. The cluster peak is not well iso- 
lated in the velocity space in its high velocity limit, where a few 
small peaks are present. We analyze the two closest groups of 
non-member galaxies, hereafter backl andback2, formed by 15 
and 11 galaxies, respectively (see Fig.lSj. Figure [16] shows the 
spatial distribution of these "background" groups. Galaxies of 
backl and back2 have a spatial distribution different from other 
non-member galaxies (at the 99.96% c.l. 95% c.l. according to 
2DKS-test). While galaxies of backl lie at the SW edges of the 
sampled field and seem to have no connection with the cluster, 
galaxies of back2 are loosely distributed in central-intermediate 
cluster regions. Moreover, back2 group is characterized by a re- 
markably small velocity dispersion cr^ ~ 150 km s"' and a dis- 
tance of Az ~ 0.02 from A520. Thus back2 might be a very 
loose group connected with the cluster, since LSS connections 
are likely found between systems separated by Az ~ 0.02 (e.g., 
Arnaud et al. 2000). 

Finally, we discuss the head tail radio galaxy (ID 184). 
According to Bliton et al ( 1 19981 1 this class of radio galaxies 
could be related with bulk motion of the intergalactic medium 
and might indicate the presence of a merger. ID 184 has the tail 
oriented opposite to the cluster center as expected in a radio- 
galaxy infalling onto the cluster (see Figs. [T] and [TTT i. Having 
high velocity v~ 60972 km s"', ID 184 is likely not a galaxy of 
the low velocity structure £, rather it might be connected to the 
infall along the LOS. 



6. Conclusions: A520 at the crossing of three LSS 
filaments 

Our findings agree with a scenario where A520 is forming at the 
crossing of three filaments of the LSS: the NE-SW one, the E-W 
one, the one along about the LOS. Clusters are expected to form 
through the accretion along three main filaments according to the 
LSS formation in the CDM scenario (Shandarin and Zeldovich 
119891 see d so beautiful images of simulated clusters Springel 
et al. |200dn) . Indeed, a few examples of clusters forming at the 



see |http://www.mpa-garching.mpg.de/galform/millermium/| 



crossing of two filaments were already observed (e.g., Arnaud et 
al.l2000l Cortese et al. "SOOT; Boschin et al. "SOOT: Braglia et al. 
I2007I I - see also Matsuda et al. (2005 ) for a protocluster at the 
crossing of three filaments. 

In this scenario the massive dark core found by M07 anal- 
ysis would coincide with the peak of the collisional component 
as shown by X-ray data only due to the particular angle of view 
of the observer. In fact, the X-ray peak likely traces the poten- 
tial well of the forming cluster, while the filament aligned with 
the LOS, projected onto the cluster center, would produce the 
peak in the 2D mass distribution. Indeed, the hypothesis of a 
LSS filament projected onto the location of peak No. 3 was al- 
ready suggested by M07 since is not in obvious contrast with 
their gravitational lensing data and X-ray data. 

Our analysis shows how powerful is the study of the internal 
cluster dynamics on the base of velocities and positions of mem- 
ber galaxies. It provides additional information which comple- 
ments X-ray and gravitational lensing analyses. Other insights 
into A520 might be recovered from the knowledge of galaxy 
properties (see e.g. Ferrari et al. 2003; Boschin et al. 2004). In 
particular, important information comes from the spectral types 
of member galaxies, since star formation could increase or, al- 
ternatively, stop during the merging phase; thus the spectral sig- 
natures of past activity are useful to determine the relevant time- 
scales (e.g., Bekki 1999; Terlevich et al. |1999l ). We are planning 
further studies of A520 in this perspective. 
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